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We present an examination of high resolution, ultraviolet spectroscopy from 
Hubble Space Telescope of the photospheric spectrum of the O-supergiant in the 
massive X-ray binary HD 226868 = Cyg X-l. We analyzed this and ground- 
based optical spectra to determine the effective temperature and gravity of the 
09.7 lab supergiant. Using non-local thermodynamic equilibrium (non-LTE), 
line blanketed, plane parallel models from the TLUSTY grid, we obtain T e ff = 
28.0 ± 2.5 kK and log g > 3.00 ± 0.25, both lower than in previous studies. The 
optical spectrum is best fit with models that have enriched He and N abundances. 
We fit the model spectral energy distribution for this temperature and gravity 
to the UV, optical, and IR fluxes to determine the angular size of and extinction 
towards the binary. The angular size then yields relations for the stellar radius 
and luminosity as a function of distance. By assuming that the supergiant rotates 
synchronously with the orbit, we can use the radius - distance relation to find 
mass estimates for both the supergiant and black hole as a function of the distance 
and the ratio of stellar to Roche radius. Fits of the orbital light curve yield an 
additional constraint that limits the solutions in the mass plane. Our results 
indicate masses of 23lgM Q for the supergiant and III3M0 for the black hole. 

Subject headings: binaries: spectroscopic — stars: early-type — stars: individual 
(HD 226868, Cyg X-l) — ultraviolet: stars — X-rays: binaries 



1. Introduction 



The massiv e X-ray binary Cygnus X-l = HD 226868 consists of an 09.7 lab primary 
( jWalbornl Il973l ) with a black hole (BH) companion. The fundamental properties of this 
system h ave been the subject of many stud ies, but they continue to be controversial. For 
example, IShaposhnikov &: Titarchukl (120071 ) determined a relationship between black hole 
mass and observed X-ray properties in the low frequency, quasi-periodic oscillation - spectral 
index plan e to derive a BH mass of 8.7±0.8M M for Cyg X-l , a value at the low end of previous 
estima tes (jGies &: Bolton! Il986al : lAbubekerov et al.l 120051 ) . On the other hand, IZiolkowski 
( 120051 ) used temperature - luminosity relations in conjunction with evolutionary models to 
calculate t h e mas s of the bright, mass donor sta r. Then, w i th the orbital mass function from 
Gies et al.l (120031 ) and the method outlined by iPaczyhskil (119741 ). he estimated the mass of 
the BH as 13.5 — 29M , at the high end of prior estimates. 

Our goal in this paper is to determine if mass estimates from these two methods can be 
reconciled through a re-examination of the supergiant 's spectrum to determine the stellar 
temperature, mass, and radius. Shortly after the X-ray source Cyg X-l was identified with 
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the star HD 226868 flBoltonl Il972l ; IWebster & Murdinl Il972f ). IWalborrJ (119731 ) classified it 
as a normal 09.7 lab star. The stellar temperature of a star of this type dep ends critically 



upon the model atmosphere assumptions adopted to match the line spectrum ([Martins et al. 



2005h. From a classi cal curve of growth analysis of the optical spectrum of HD 226868, 
Canalizo et al.l (119951 ) estimated an effective temp erature of T e w = 32 ± 2 kK, and they 



found an overabundance of He in the photosphere. iHerrero et al.l (119951 ) also estimated the 
temperature of the star T eff ps 32 kK based upon fits of the optical line spectrum with cal- 
culated profiles from unified model atmospheres that included a non-LTE treatment of H 
and He but neglected line-blanketing from transitions of heavier elements. In addition, they 
determined values for gravity from fits of the Balmer lines that ranged from log g = 3.03 for 
plane-parallel models to log g = 3.21 for spherical models that included wind effects. Their 
results led to mass estimates of 17.8 a nd 1O.1M for the supergiant and BH, respectively. 
More recently, iKaritskaya et al.l (120081 ) classified HD 226868 as an ON star with a temper- 
ature of T e fj = 30.4 ± 0.5 kK and gravity of \ogg = 3.31 ± 0.07 using a semi-gray model 
atmosphere that accounts for non-LTE effects in some lines and for X-ray illumination. 

Here we present an analysis of the photospheric parameters for the supergiant based 
upon ground-based optical spectra and high-resolution, UV spectra from the Hubble Space 
Telescope Sp_ace Telescope I magi ng Spectrograph (S TIS). These STIS spectra were first pre- 
sented by iGies et al.l (120081 ) and IVrtilek et al.l (120081 ) in discussions of the orbital variations 
observed in the stellar wind lines. We compare the optical and UV line profiles of HD 226868 
with synthetic spectra based on line blanketed, non-LTE photospheric models in order to 
determine the stellar temperature and gravity (§2). Since the continuum flux and spectral 
lines of the supergiant could be influenced by X-ray heating, we search for heating effects 
in the orbital UV flux variations using the low/hard state International Ultraviolet Explorer 
(IUE) archival spectra and the high/soft state HST spectra (§3). A stellar radius - distance 
relation can be determined from fits of the spectral energy distribution. We use the observed 
flux distribution and spectra of field stars in the same region of the sky to estimate the 
reddening and extinction in the direction of Cyg X-l and to determine the angular size of 
the star (§ 3). Fi nally, we use this radius - distance relation with the method developed by 



Paczyhskil (119741 ) to set mass limits as a function of distance and to estimate the probable 
masses using constraints from the rotational line broadening and ellipsoidal light curve (§4). 



2. Ultraviolet and Optical Line Spectrum 



We need to rely on the line spectral features to estimate temperature since the UV and 
optical continuum falls in the long wavelength, Rayleigh- Jeans part of the flux distribution 
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where the shape of the continuum is insensitive to temperature. Some of the best line 
diagnostics for late O-supergiants are found in the optical spectrum where several ionization 
state lin e ratios and the Balmer line profiles change dr amatically with temperature and 



gravity (I Walborn fe Fitzpatrickl Il990l ; ISearle et al.l 120081 ) . In this section, we use xl fits 



of the optical and UV spectra with model spectra to estimate T e s and logg. Our data 
consist of high resolution UV spectra taken with the HST/STIS (G140M grating, resolving 
power R = 14500) and two sets of optical spectra from the Kitt Peak National Observatory 
Coude Feed telescope (CF; 3759 - 5086 A, R = 2990) and 4 m Mayall Telescope and RC 
spectrogra ph (RC; 4182 - 49 42 A, R = 5700). Details of these observations are given in 
Table 1 of iGies et al.l (120081 ). All of thes e flux-rectified sp ectra were shifted to the rest 
frame (using the orbital solution given by IGies et al.l 120081 ) and co-added to increase the 
signal-to-noise ratio. 

We compared these spectra with model spectra fr om the TLUSTY/ SYNSPEC codes 
given in the grids OSTAR2002 JLanz fe Hubenvi l2003h and BSTAR2006 (jLanz fe Hubenv 
20071 ). The model atmospheres are based upon a plane parallel geometry, solar abundances, 
line blanketed opacities, and non-LTE calculations of atomic populations for H, He, and 
representative atoms up to Fe. The model spectra are presented as a function of four param- 
eters: the microturbulent velocity of gas in the line forming region, £, the stellar effective 
temperature, T e g, logarithm of the gravitational acceleration in the photosphere, log g, and 
the chemical abundance of the gas. The model spectra were transformed to the observed 
wavelength grids by wavelength integration and convolution with rotational and instrumen- 
tal broadening function s. We adopted a projected rotational velo city of V sin i = 98 km s" 1 



( IGies fe Boltonl Il986al ). linear limb darkening coefficients from IWade fe Rucinskil (119851 ). 
and Gaussian rep resentations of instrumental broadening using the projected slit FWHM 
jGies et alJliooi ). 



It was clear from inspection that a large microturbulence is required to match the 
observed and model spectra. The OSTAR2002 grid uses £ = 10 km s _1 throughout while 
the BSTAR2006 grid adopts £ = 2 km s _1 for the full grid and £ = 10 km s -1 for a selection 
of low gravity (supergiant) models. We found that the best fits were obtained in all three 
spectral bands with the £ = 10 km s _1 models, and this was especially true in the FUV where 
the observed deep spectral lines were not matched with the lower microturbulent velocity 
models. An atmosp heric microturbulence of £ = 10 km s _1 is typical for late-0 supergiants 
( iRyans et al.ll2002l ). and lCanalizo et al.l (119951 ) derived an estimate of £ = 10.7 km s _1 from 
a curve of growth analysis of the optical N III lines in the spectrum of HD 226868. 



We tested the goodness-of-fit for each of the models of interest by calculating the reduced 
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xl statistic, 



A? 




(i) 



i=0 



Here N is the total number of wavelength points used in the fit and <r err (Aj) is the standard 
deviation of the mean rectified flux (determined from the scatter at wavelength Aj among 
the individual spectra in the co- added mean). We selectively omitted from the summation 
spectral regions that contained stellar wind features or interstellar lines that are not present 
in the model spectra. Our results are listed in Table 1 for a wide range of model spectra with 
a microturbulence of £ = 10 km s _1 . Column (1) indicates the spectral region fit by "HST" 
for the FUV spectrum, "CF" for the KPNO CF, blue spectrum, and "RC" for the KPNO 
4 m RC, green spectrum. Column (2) gives the grid value of gravity logg, and column (3) 
gives a code for the spectral model ("O" for spectra from the OSTAR2002 grid and "B" 
and "BCN" for spectra from the BSTAR2006 grid). Then follow 10 columns that list the 
measured xl f° r grid values of T cff (at increments of 2.5 kK and 1 kK for the OSTAR2002 
and BSTAR2006 grids, respectively). 

The trends in Table 1 are represented in a combined contour diagram in Figure 1. Here 
the gray-scale contours represent the goodness-of fit for the FUV spectrum, the solid lines 
for the blue spectrum, and the dashed lines for the green spectrum. Since there is not an 
exact match between the predictions of the OSTAR2002 and BSTAR2006 grids at their 
boundary, Figure 1 shows contours based only on the OSTAR2002 grid for high gravity 
models logg > 3.0, while the contours in the lower gravity region log g < 3.0 are based upon 
the BSTAR2006 grid. Note that there are no models available for high temperature, low 
gravity region in the lower right part of the diagram because such atmospheres approach 
or exceed the Eddington luminosity limit. The higher gravity, lower temperature region 
is empty because the BSTAR2006 grid contains only models for a lower microturbulent 
velocity £ = 2 km s _1 for this parameter range. Note that the xl minima in Table 1 have 
values much larger than the expected value of unity. This is due to the inclusion of spectral 
regions where there is evident mismatch because of incomplete removal interstellar features, 
marginal differences in the continuum placement, and real differences between the observed 
and models even in the best fit cases. For the purposes of intercomparison in Figure 1, we 
have subtracted from each xl set the minimum minus one value, so that the figure represents 
variance increases that result only from changes in the assumed temperature and gravity. 

The temperature and gravity properties of the xl fits shown in Figure 1 result primarily 
from the dependence of the spectral lines on the ionization levels in the gas. The Saha ioniza- 
tion equilibrium equation shows that the number ratio of atoms of one ion to the next higher 
ionization state is equal to the electron density times a function that decreases with increas- 
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ing temperature. Thus, in order to match the ionization ratios represented in the spectral 
line strengths, the best fits are found along a diagonal zone of increasing log g with increasing 
temperature, i.e., increasing the electron density with logg compensates for the functional 
drop related to temperature increase. However, there are other spectral dependences that 
help us determine the minimum xl position along the valley in the (T e s,\ogg) diagram. In 
particular, the H Balmer line wings in the optical spectrum are sensitive to pressure broad- 
ening (linear Stark effect) and hence model gravity. The lowest contours of xl i n Figure 1 
indicate that the best fits are found for T eff = 26.5 — 28.5 kK and log.q = 2.9 — 3.1. Both of 



these estimates are lower th an found in earlier studies (jCanalizo et al.l Il995l ; iHerrero et al. 



19951 ; iKaritskaya et al.l 120081 ) . but they are consistent with recent studies that demonstrate 
that the inclusion of li ne blanketing in stellar atmosphere models tends to lower the derived 



effect ive temperature (IRepolust et al.l 12004 ; iMartins et al.l 120051 ; iLefeverl 120071 ; ISearle et al. 
2008h . 



The run of xt fi^ s shown in Figure 1 are for solar abundance models, but given the 
reports that the spectrum of HD 226868 has strong N lines, we also explored fits with 
CN abundance models included in the BSTAR2006 grid for gravities of log g = 2.75 and 
3.00. These models assume a He to H number ratio of 0.2 (compared to 0.085 for the solar 
abundance models), a C abundance equal to one half the solar value, and a N abundance 
five times the solar abundance. These adjustments demonstrate the kind of changes that 
are expected when the atmosphere becomes enriched in CNO-processed gas. While the CN 
models do not improve the fits of the FUV spectrum, they are significantly better fits of 
the optical spectrum (where a number of strong He I and N III lines are present; see Fig. 2 
below). When we compare the xt fit of a solar model to a CN model at the same gravity 
in Table (U the CN models fit better at higher temperature, especially for fits of the optical 
spectra. Following this trend, we estimate that the supergiant's spectrum is fit best by the 
CN models with T e g = 28 ± 2.5 kK and log g = 3.00 ± 0.25 dex. Thus, we will focus on these 
CN models for the rest of this section. 

We compare the two mean optical spectra (CF and RC) with our best fit model spectrum 
in Figure 2 (T efT = 28 kK, log g = 3.00) and with a marginally acceptable fit in Figure 3 (T cff = 
26 kK, log g = 2.75). Als o shown for comparison is the spectrum of a similar 09.7 lab star, 
\i Nor (HD 149038; from lWalborn fc Fitzpatrickj|l990l ). The spectra appear in three plots: 
the short wavelength range from the lower resolution CF data is shown in the top panel while 
the bottom two panels illustrate the longer wavelength region from the higher resolution RC 
spectra (compare with Fig. 1 in IKaritskaya et al.ll2008l ). In general the agreement between 
the observed and model spectrum is satisfactory. The largest discrepancies are seen in the 
He II A4686 and B.0 lines where incipient emission from the stellar wind of t he su pergiant 
alters the profiles Jdes fc Boltonill986bl ; iNinkov et aDll987al : bies et al.lfiooil . l2008h . The H 
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Balmer line emission is strongest in Ha, and for simp le estimates of the Balmer decrement for 
Case B recombination flOsterbrock &: Ferlandll2006l ). we expect some measurable degree of 
wind emission for all the Balmer lines shown in Figures 2 and 3. We find that the H line cores 



do ap pear shallower, while the Balmer line wings agree well with the models. ISearle et al. 



(120081 ) observed this effect in other supergiants and they suggest that stellar wind emission 
from the outer atmosphere tends to fill in the line core. On the other hand, the Balmer 
line wings are formed in higher density gas, deeper in the photosphere where we expect the 
TLUSTY/SYNSPEC results to be quite reliable. The H line wings become narrower with 
lower gravity, and the predicted H profiles for the lower gravity model illustrated in Figure 3 
appear to be significantly narrower than the observed ones. 

The He I and He II line strengths are well matched in the He enriched CN model 
spectra. In particular the temperature sensitive ratio of He II A4541 to Si III A4552 
(equal for the 09.7 classification) is better reproduced by the T e g = 28 kK and log g = 
3.00 model (Fig. 2) than the T cS = 26 kK and \ogg = 2.75 model (Fig. 3). The N III 
AA4097, 4379, 4510, 4514, 4630, 4634, 4640 lines are also reasonably well fit in the five times 
overabundant CN models. On the other hand, O lines like O II AA4069, 4072, 4075, 4590, 4596 
are too strong in the model spectra, which suggests that the O abundance should be revised 
downwards from solar values as expected for CNO-processed gas. The other differences be- 
tween the observed and model spectra are related to the presence of interstellar features 
(Ca II AA3933, 3968 in the CF spectrum; most of the deep ISM features were removed from 
the RC spectra). 

In Figure 4 we present the averaged UV spectrum made with HST/ST1S with the best 
TLUSTY/SYNSPEC model superimposed as a lighter line. Figure 4 also includes an average 
UV spectrum of /i Nor, based upon 34 high resolution, archival spectra from IUE. Horizontal 
line segments indicate those regions where the lines primarily originate in the photosphere, 
i.e., free from P Cygni st ellar wind lines and from regions where interstellar lines were 
removed by interpolation (IGies et al.ll2008l ). Overall, the line features in the observed UV 
spectrum agree well with the model UV spectrum based upon the optimal T e s and log g 
parameters derived from the optical and FUV spectral fits. Note that the He II A 1640 
feature appears in absorption as predicted, s o the re is no evidence of the Raman scattering 
emission that was observed by iKaper et al.l (Il990l ) in the massive X-ray binary 4U1700-37. 
There are, however, a few specific regions where the match is less satisfactory. For example, 
the blends surrounding Fe V A 1422 and Fe IV AA1596, 1615 appear stronger in both the 
spectra of HD 226868 and \x Nor, which suggests that the mo dels are underes timating the Fe 
line opacity in these wavelength regions. The S V A 1502 line (jHowarthl 119871 ) has a strength 
in the spectrum of HD 226868 that falls between that of the model and of \x Nor. The deep 
feature near 1690 A is an instrumental flaw near the edge of the detector at one grating tilt. 
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There are huge variations in the stellar wind lines between the orbital con 



unctions 

that are due to X-ray ionization of the wind (IGies et al.l 120081 ; IVrtilek et al.ll2008l ). and it 
is possible that X-ray heating might also affect some of the photospheric lines. Figure 5 
compares the average UV spectra at the two conjunction phases <p — 0.0 and 0.5 (inferior 
and superior conjunction of the supergiant, respectively). With the exception of the known 
wind line changes, we find that the spectra are almost identical between conjunctions. Some 
slight differences are seen in very strong features, such as the Si III A1300 complex and the 
Fe V line blends in the 1600 - 1650 A region. The deep lines appear somewhat deeper 
at cf) = 0.0 and have slightly extended blue wings compared to those observed at <fi = 0.5 
(when the black hole is in the foreground). We speculate that the deeper cores and blue 
extensions result from line opacity that forms in the upper atmosphere where the outward 
wind acceleration begins. This outer part of the atmosphere in the hemisphere facing the 
black hole may also experience X-ray ionization (like the lower density wind) that promotes 
Si and Fe to higher ionization levels and reduces the line opacity of the observed transitions. 

Our spectral fits are all based upon the existing OSTAR2002 and BSTAR2006 grids, 
and it would certainly be worthwhile to explore more specific models, for example, to derive 
reliable estimates of the He and N overabundances. A determination of the He abundance 
in particular will be important for a definitive temperature estimate. It is als o important in 



such a n analysis to consider the full effects of the stellar wind in HD 226868. iHerrero et al. 



( 119951 ) compared analyses of the spectrum of HD 226868 from static, plane-parallel models 
with unified, spherical models (that treat the photosphere and wind together), and they 
found their log g estimate increased by about 0.2 dex (with no change in temperature) in 
the unified models. Thus, we suspect that our gravity estimate derived from the plane- 
p arallel TLUSTY cod e is probably a lower limit (approximately consistent with the results 



of Herrero et al. 1995 and 



Karitskaya et al.ll2005l ). 



3. UV - IR Spectral Energy Distribution 



We can use the derived model flux spectrum to fit the observed spectral energy dis- 
tribution (SED) and reassess the interstellar extinction and the radius - distance relation. 
We collected the archival low dispersion IUE spectra and combined these fluxes with the 
HST spectra in wavele ngth bins spanning t he FUV and NUV regions. We tra nsformed the 
UBV magnitudes from iMassey et al.l (119951 ) into fluxes using the calibration of IColina et al. 
( 119961 ). and the near-IR fluxes were determined from a calibration of the 2MASS J UK, 



magnitudes (jCohen et al.ll2003l ; ISkrutskie et al.ll2006l ). Then we fit the observed fluxes with 
the optimal BSTAR2006 flux model (CN model, f = 10 km s" 1 , T eff = 28 kK, log# = 3.0) to 
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find the best reddening curve using the extinction law from iFitzpatrickl (119991 ) . We placed 
additional weight on the six optical and IR points to compensate for the larger number of 
UV points. Figure 6 shows the observed and best fit model fluxes for HD 226868 that we 
obtained with a reddening E(B — V) — 1.11 ± 0.03 mag and a ratio of total to selective 
extinction Ry = 3.02 ± 0.03. These values agree well with the previous reddening estimates 
that are collected in Table 2. 

For comparison we examined the colors and reddening of six field stars within 10 ar- 
cminutes of HD 226868 in the sky. These stars were observed with the KPNO 4 m telescope 
and RC spectr ograph using the same blue region arrangement selected for our observations 
of HD 226868 (IGies et al.ll2008l ). We made spect r al cla ssification of the stars, and then we 
used the observed UBV colors fro mlMassev et al.l (119951 ) and the intrinsic color and absolute 
magnitude for the classification (lGrayill992l ) to estimate reddening and distances to these 
stars. Our r esults are collec t ed in Table 3 with the reddening estimate from above listed for 
HD 226868. bregman et"aD Jl973h estimated the distance of HD 226868 as d « 2.5 kpc, and 
set a lower limit of 1 kpc based upon the colors of other nearby field stars. We find that 
there are two stars at distances just under 1 kpc that have a similar reddening to that of 
HD 226868, which is consistent with a distance to HD 226868 of d > 1.0 kpc. 

The normalization of the fit to the SED yields the star's limb-darkened angular diameter, 
9 = 96 ± 6 yuas. Then we can calculate the luminosity and radius of the star as a function of 
distance d (in kpc) to HD 226868, 

h. 



(5.9 ±2.1) x 10 4 d 2 , 



R 

- 1 = (10.3 ± 0.7)d. 

R e 



(2) 
(3) 



It is important to check that these SED results are not affected b y long term or or bital 
flux variability, so we examined the archival IUE low dispersion spectra (IGies et al.ll2008l ) and 
the HST spectra to determine the amplitude of any flux variations in the UV. We calculated 
the average continuum flux over three wavelength spans (1252 - 1380 A, 1410 - 1350 A, and 
1565 - 1685 A) that excluded the main wind features. We then converted the UV fluxes 
to differential magnitudes Am. We found no significant differences between fluxes from 
tim es corresponding to the X-ray low/hard state (IUE) and high/soft state (IUE and HST] 



sec 



Gies et al.l 120081 for X-ray state information), nor were there any long term variations 
over the 25 year time span between the IUE and HST observations. On the other hand, 
we do find marginal evidence of the orbital flux variations related to the tidal distortion of 
the supergiant. We plot in Figure 7 the mean orbital flux variations of the three wavelength 
intervals for both IUE and HST spectra that are averaged into eight bins of orbital phase. For 
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compa rison we also include the \^-band ellipsoidal light curve from iKhaliullin fc Khaliullina 
()198ll ). The UV an d l/-band light curves appear to have s imilar amplitudes, consistent 
with past estimates ([Treves et al.lll980l ; Ivan Loon et al.ll200ll ). N ote that the minim a have 
approximately equal depths (consistent with the optical results; iBalog et al.lll98ll ). which 
suggests that there is little if any deep heating by X-rays of the hemisphere of the supergiant 
facing the black hole. Since the amplitude of the light curve is small and the average UV 
fluxes plotted in the SED in Figure 6 cover the full orbit, the ellipsoidal variations have a 
minimal impact on the quantities derived from the SED. 

Finally, we need to consider if the SED has a non-stellar flux contribution from the 
accretion disk around the black hole or from other circumstellar gas. iBruevich et al.l (119781 ) 
estimated that the disk contributes about 2% of the optical flux, and there are reports of 
small optical variatio ns with superorbital periods that may correspond to the precession 
of the accretion disk (IKemp et al.l 119871; iBrocksopp et al.l Il999l ; ISzostek fc Zdziarskil 120071 ; 
Poutanen et al.ll2008l ). Furthermore, iDolanl (|200ll ) observed rapid UV variations that he 
argu ed originate i n dyin g pulse trains of infalling material passing the event horizon of the 
BH. [Miller et al.l (120021 ) developed a multi-color disk SED to model the X-ray continuum 
of Cyg X-l, and Dr. Miller kindly sent us the model fluxes extrapolated into the UV and 
optical. These are also plotted in Figure 6 after accounting for interstellar extinction. Both 
the photospheric and disk SEDs correspond to the Rayleigh- Jeans tail of a hot continuum, 
and the model predicts that the disk contributes approximately 0.01% of the total flux in the 
UV to IR range. This small fraction is consistent with our successful fitting of the UV and 
optical line features that would otherwise appear shallower by flux dilution if the disk was a 
significant flux contributor. Thus, our SED fitting is probably unaffected by any non-stellar 
flux source. 



4. Mass of the Supergiant 



In this section we will explore the mass c onseq uences of our relations for radius and 
luminosity as a function of distance. IPaczyhskil ( 1974 ) derived model- independent, minimum 
mass estimates for both components as a function of distance based on the lack of X-ray 
eclipses (setting a maximum orbital inclination) and the assumption that HD 226868 is not 
larger than its Roche lobe (setting a lower limit on the ratio of the supergiant to black hole 
mass, Mi/M 2 ). We repeated his analysis using our revised radius - distance relationship 
(eq. [3]), stellar effective temperature T c s = 28 kK, and cur rent values for th e mass function 
f(m) = 0.251±0.007 M and period P = 5.599829 days JGies et al.l liooih . The resulting 



minimum masses are presented in columns 7 and 10 of Table 4 as a function of distance d. 
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We can make further progress by assuming the supergiant has attained synchronous 
rotatio n with the orbit since the stellar radius is probably comparable in size to the Roche 
radius (IGies fc Bolton! Il986bl ). We take the ratio Q of the star's spin angular velocity to 
orbital angular velocity to be 1. Then the projected rotational velocity Vsini is related to 
the inclination i by 



V sin i 



2tt 
~P 



Ri sin i 



(4) 



where P is the orbital period. The projected rotational velocity, afte r correction for macro- 

r sini = 95 ± 6 km s" 1 jGies fc Bolton] Il98( 
iserting equation Q for R\ we obtain an inch 



turbulent broadening, is estimatec 


1 to b 


Ninkov et al. 


1987b 




Herrero et al. 


1995) 



nation estimate in terms of distance d (kpc) of 

i = arcsin((1.02 ±0.09)d). 



(5) 



These inclination estimates are given in column 2 of Table 4. Note that this argument 
suggests a lower limit to the distance of ~ 1.0 kpc, similar to that found by reddening 
considerations. 



Gies &: Bolton! (11986a! ) showed that the mass ratio can be estimated from the ratio of 
the projected rotational velocity to the orbital semiamplitude of the supergiant, 



V^sini 
K 



p(Q + l)$(Q) 



(6) 



where p is the fill-out factor, i.e., the ratio of volume equivalent radii of the star and Roche 
lobe, Q = My /My, is the mass ratio, $ is the ratio of the Roche lobe radius to the semimajor 
axis (!Eggletonlll983l ). and synchronous rotation is assumed. Thus, given the observed values 
of Vsini and K and an assumed value of p, we can find the mass ratio and, with the 
inclination, the masses of each star. These masses are listed in columns 8, 9, 11 and 12 of 
Table 4 under headings that give the fill-out factor in parentheses. The run of masses is also 
shown in Figure 8 that illustrates the mass solutions as a function of distance and fill-out 
factor. Loci of constant p are denoted by dotted lines (increasing right to left from p = 0.85 
to 1.0) while loci of constant distance (and inclination angle) are shown by dashed lines. The 
derived gravity values from these masses of logg ~ 3.3 reinforces the idea that our spectral 
estimate of logg = 3.0 is a lower limit (see §2). 

We assumed synchronous rotation in the relations above because both observations and 
theory indicate that the orbit al synchronization time scale in close binaries is shorter than 
the circularization time scale (IClaret et al.lll995l ). and since the orbit is circular, it follows 
that the star must rotate at close to the synchronous rate. However, it is straight forward 
to see how the solutions will change if the synchronism parameter Q differs from unity. In 
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equation (5) the distance d can be replaced by the product Qd, while in equation (6) the 
fill-out parameter p can be replaced by Qp. If, for example, Q = 0.95, then the mass solutions 
can be obtained from Table 4 and Figure 8 by selecting a distance of 0.95<i and a fill-out 
ratio of 0.95p. 

The other important constraint comes from the ellipsoidal light curve. The tidal distor- 
tion of the star results in a double-wave variation (Fig. 7) whose amplitude depends on the 
inclination (maximal at % = 90°) and degree of tidal distortion (maximal for fill-out p = 1.0). 
In order to determine which parts of mass plane are consistent with the observed variation, 
we constructed model IZ-ba nd light curves using the GENSYN code (IMochnacki fc Doughty 
19721 ; iGies fc Boltonlll986al ) for the four values of fill-out factor illustrated in Figure 8. There 
is a unique solution for the best fit of the light curve along each line of constant fill-out 
factor, since the light curve amplitude monotonically decreases with decreasing inclination 
(increasing distance). The solid line in Figure 8 connects these best fit solutions (indicated 
by plus sign symbols) . These light curve solutions differ slightly from those presented by 
Gies fc Bolton! (1l986aj ) beca use we chose t o fit t he light curve from iKhaliullin fc Khaliullina 
(1l98ll ) instead of that from iKemp et al.l (119831 ). and the differences in the solutions reflect 
the uncertainties in the observed light curve. 

There are several other constraints from hints about the mass transfer process, lumi- 
nosity, and distance th at c an provide additional limits on the acceptable mass ranges. Both 
Gies &: Bolton! (Il986bl ) and lNinkov et al.l (ll987al ) presented arguments that the unusual He II 
A4686 emission in the spectrum of HD 226868 originates in a tidal stream o r focuse d wind 
from the supergiant towards the black hole. Furthermore, IGies fc Bolton! (jl986bl ) made 
radiative transfer calculations of the focu sed wind emission profiles for models of the asym- 
metric wind from lFriend fc Castor! (119821 ). and they determined that the fill-out factor must 
exceed p = 0.90 in order to increase sufficiently the wind density between the stars to account 
for the observed strength of the He II A4686 emission. Thus, the presence of a focused wind 
implies that the fill-out factor falls in the range p = 0.9 — 1.0. 



Paczyhskil (11974] ) and IZiolkowskil (120051 ) argue that massive stars evolve at near con- 



stant luminosity, and, therefore, the best solutions will obey the observed mass - luminosity 
relation. Table 4 lists the derived luminosity log Li (column 4) as a function of distance (eq. 
[2]) plus the predicted luminosities for the mass solutions determined for the p = 0.9 and 
1.0 cases, logL*(0.9) and logL^(l.O), respectively (columns 5 and 6). These predictions are 
based upon the mass - lu minosity relations for T e s = 28 kK stars from the model evolu- 
tionary sequences made by ISchaller et al.l (jl992j). We find that the observed and predicted 
luminosities match over the distance range of d = 1.7 (p = 0.9) to 2.0 kpc (p = 1.0), closer 
than the range advocated by lZiolkowskil (120051 ) who adopted a higher temperature and hence 
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higher luminosity. Note that some stars in mass transfer binaries appear overluminous for 
their mass, so these distances should probably be considered as upper limits. 

Several authors have suggested that the position an d proper motion of HP 2 26868 
indicates that it is a member o f the Cyg OB3 associ ation (IMirabel fc Rodriguesl 120031 ) that 
has a distance of 1.6 -2.5 kpc (jUyaniker et al.l 120011 ). However, a radio parallax study by 
Lestrade et al.l (1l999l ) indicates a smaller (but possibly consistent) distance of 1.4^4 kpc 
for Cyg X-l. Our fits of the ellipsoidal light curve suggest that the maximum allowable 
distance is d ~ 2.0 kpc (for p = 1.0) The interstellar reddening indicates a distance of at 
least 1.0 kpc (§3), w hich is probably consistent with the strength of interstellar Ca II lines. 
Megier et al.l (120051 ) present a method for determining the distance to O supergiants using 
the equivalent width W\ of t he Ca II A3933 f eature . Using their calibration with the value 
of W\ = 400 ± 10 niA from iGies fc Bolton! (11986a! ) yields a distance d = 1.2 kpc. Since 
the reddening of H P 226868 is approxim ately the same as that for the much more distant 
Cepheid, V547 Cyg (IBregman et al.lll9731 ). the ISM must have a relatively low density beyond 
m 1 kpc along this line of sight through the Galaxy, so we suspect that the distance derived 
from the interstellar Ca II line is probably a lower limit. 

All of these constraints are consistent with the mass solutions for a fill-out factor range 
of p = 0.9 — 1.0, and the corresponding mass ranges are listed in Table 5. We also list mass 
estimates from earlier investigations. Our downwar d revision o f the e ffective temperature re- 
sults in lower luminosity estimates than adopted by lZiolkowskil (120051 ). and consequently, our 
mass estimates (based upon the light curve) are significantly lower than his mass estimates 
(based upon the mass - luminosity relation from models). In fact, the lower limit for the black 
hole m ass now overlaps comfortably with the mass determined by lShaposhnikov fc Titarchuk 
(120071 ) using the correlation between the X-ray quasi-periodic oscillation frequency and spec- 
tral index, so the apparent discrepancy in black hole mass estimates from X-ray and optical 
data is now resolved. If the X-ray derived mass is accurate, then the mass solution for fill-out 
factor p = 0.91 is preferred (Mi = 19M and the distance is d = 1.6 kpc). 

Our analysis of the first high resolution UV spectra of HP 226868 and of the com- 
plementary optical spectra shows that the photospheric line spectrum can be matched 
by adopting an atmosphere mixed with CNO-processed gas with an effective temperature 
T e fj = 28.0 ± 2.5 kK and log g > 3.0 ± 0.25. Assuming synchronous rotation (Q = 1) 
and using the fill-out factor range from above, the mass of the supergiant ranges from 
Mi = 17 — 31M and the black hole mass ranges from M 2 = 8 — 16M Q . This corresponds 
to an inclination of i = 31° — 43° and a distance of d = 1.5 — 2.0 kpc. Better estimat es of 
the masses may be possible i n the future. For example, both the GAIA (I Jordan! 120081 ) and 
SIM Lite ( Unwin et ali 2008 ) space astrometry missions will provide an accurate parallax 
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and distance. Furthermore, pointed observations with SIM Lite will measure the astrometric 
motion of the supergiant around the system center of mass, yielding independent estimates 
of both the orbit al inclination and dist ance (by equating the astrometric and radial velocity 



semi-major axes; iTomsick et al.l 120091 ) . Finally, future high dispersion X-ray spectroscopy 
with the International X-ray ObservatorS^ will measure the orbital motion of the black hole 
through the orbital Doppler shifts of accretion disk flux in the Fe Ka line (IMillerl 120071 ). By 
comparing the optical and X-ray orbital velocity curves, we will have a secure mass ratio that, 
together with the distance estimate, will lead to unique and accurate mass determinations 
of the supergiant and black hole. 
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Table I. xl for Spectral Fits with Models 
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Table 2. Interstellar Reddening Estimates 



Source 


E(B 


-V) 


Ry 


Brctrman et al. (1973) 


1.12 


(5) 


3.0 


Treves et al. (19801 


1.06 




3.0 


Wu et al. (1982) 


0.95 


(7) 


3.1 


Savage et al. (1985) 


1.080 (25) 


3.1 


Weener (2002) 


1.03 




3.32 


This Paper 
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Note. — Numbers in parentheses give the error 
in the last digit quoted. 



Table 3. Interstellar Reddening for HD 226868 and Nearby Stars 
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Table 4. Mass and Luminosity versus Distance for HD 226868 
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Table 5. Mass Estimates 





Mi 


M 2 


Source 


(M @ ) 


(M ) 


Baloe et al. (19811 


20 - 27 


7-12 


Gies & Bolton fl986a1 


23 - 38 


10 - 20 


Ninkov et al. f 1987b 1 


20 


10 


Herrero et al. f 19951 


17.8 


10.1 


Abubekerov et al. (20051 


22 


8.2 - 12.8 


Ziolkowski f 20051 


30 - 50 


13.5 - 29 


Shaposhnikov & Titarchuk f2007l 




7.9 - 9.5 


This paper 


17 - 31 


8-16 
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Fig. 1. — The variation within the (T e g,\ogg) plane of the net reduced xl statistic that 
measures the goodness of fit between the solar abundance models and the observed spectrum 
of HD 226868. The contours show the value of xt above the best fit minimum (arbitrarily 
set to one), and they nominally represent the intervals of 2a, 4a, ...10a where a is the error in 
the parameter estimate. The contours for the FUV, KPNO CF, and KPNO 4 m spectral fits 
are shown as different gray-shaded regions, solid lines, and dashed lines, respectively. The 
contours for log g < 3.0 are for fits with the BSTAR2006 models, while those for \ogg > 3.0 
are for fits with the OSTAR2002 models. 
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Fig. 2. — Rectified optical spectra (dark line) together with a T eff = 28 kK and logg = 3.0 
TLUSTY CN model (light line). The top panel shows the mean spectrum obtained with 
the KPNO Coude Feed telescope while the bottom two panels show the mean spectrum 
obtained with the KPNO 4-m telescope. The spectrum of the 09.7 lab star /i Nor is offset 
by 0.15 from the model and HD 226868 spectra for comparison. The horizontal lines below 
the spectra indicate the wavelength regions included in the xl calculation. The He II A4686 
emission line originates in the focused wind from the star. 
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Fig. 3.— Optical spectra (dark line) together with a T eff = 26 kK and hgg = 2.75 TLUSTY 
CN model (light line) in the same format as Fig. 2. 
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Fig. 4. — Observed UV spectrum (dark line) plus a TLUSTY CN model spectrum for 
T eff = 28 kK and logg = 3.0 (light line). The spectrum of the 09.7 lab star \x Nor is offset 
by 0.6 from the model and HD 226868 spectra for comparison. The horizontal lines below 
the spectra indicate regions included in the xl calculation. 
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Fig. 5. — Averaged UV spectra of HD 226868 at orbital phases = 0.0 (light line) and 
4> = 0.5 (dark line). The horizontal line indicates those regions without strong ISM features. 
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Fig. 6. — The spectral energy distribution of HD 226868 (plus signs) together and the 
TLUSTY best fit (solid line) for T c g = 28 kK and log g = 3.0. The UV points were binned 
from the average HST and IUE spectra. The three optical points are the UBV measurements 
from iMassey et al.l (119951 ) and the three IR point s are taken from 2M ASS. Also shown is the 
extrapolation of the accretion disk flux model of Miller et al.l (120021 ) (dotted line). 
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Fig. 7. — UV light curve from both t he HST and IUE spectra {diam onds] compared with 
the l^-band light curve {circles) from iKhaliullin fc Khaliullinal (jl98ll ). The UV data were 
divided into eight orbital phase bins, and the error bars indicate the standard deviation 
within each bin. 
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Fig. 8. — Mass plane diagram for the optical (x-axis) and X-ray (y-axis) components assum- 
ing synchronous rotation of the supergiant. Dotted lines represent fill-out factors (p) of 0.85 
to 1.0 in increments of 0.05 from right-to-left, and the dashed lines show loci of constant dis- 
tance (in kpc as labeled on the left side while the corresponding orbital inclination rounded 
to the nearest degree appears on the right side). Plus si gns connected by a solid line s how 
the solutions that match the y-band orbital light curve (IKhaliullin fc Khaliullinalll98ll ) at 
each value of the fill-out factor. 



